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Abstract 

The object of this paper is the uniqueness for a d-dimensional Fokker-Planck type 
equation with inhomogeneous (possibly degenerated) measurable not necessarily 
bounded coefficients. We provide an application to the probabilistic representation 
of the so-called Barenblatt's solution of the fast diffusion equation which is the par- 
tial differential equation dtu — di x u m with m G]0, 1[. Together with the mentioned 
Fokker-Planck equation, we make use of small time density estimates uniformly with 
respect to the initial condition. 
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1 Introduction 

The present paper is divided into three parts. 

i) A uniqueness result on a Fokker-Planck type equation with measurable non-negative 

(possibly degenerated) multidimensional unbounded coefficients. 

ii) An application to the probabilistic representation of a fast diffusion equation. 

hi) Some small time density estimates uniformly with respect to the initial condition. 

In the whole paper T > will stand for a fixed final time. In a one dimension space, 
the Fokker-Planck equation is of the type 

f d t u{t,x) = d* x (a(t,x)u(t,x))-d x (b(t,x)u(t,x)), t e]0,T], x e R, 
1 "(0,-) = V{dx), 
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where a, b : [0, T] x R — > R are measurable locally bounded coefficients and ^ is a 
finite real Borel measure. The Fokker-Planck equation for measures is a widely studied 
subject in the literature whether in finite or infinite dimension. Recent work in the case 
of time-dependent coefficients with some minimal regularity was done by (9j 1161 13011 in 
the case d > 1. In infinite dimension some interesting work was produced by [8|. 

In this paper we concentrate on the case of measurable (possibly) degenerate co- 
efficients. Our interest is devoted to the irregularity of the diffusion coefficient, so 
we will set b = 0. A first result in that direction was produced in [7] where a was 
bounded, possibly degenerated, and the difference of two solutions was supposed to 
be in L 2 ([n,T] x R), for every n > (ASSUMPTION (A)). This result was applied to 
study the probabilistic representation of a porous media type equation with irregular 
coefficients. We will later come back to this point. We remark that it is not possible 
to obtain uniqueness without ASSUMPTION (A). In particular [7, Remark 3.11] pro- 
vides two measure-valued solutions when a is time-homogeneous, continuous, with - 
integrable in a neighborhood of zero. 

One natural question is about what happens when a is not bounded and x € R d . 
A partial answer to this question is given in Theorem |3.1| which is probably the most 
important result of the paper; it is a generalization of [7 , Theorem 3.8] where the in- 
homogeneous function a was bounded. Theorem |3 . 1 1 handles the multidimensional case 
and it allows a to be unbounded. 

An application of Theorem |3 . 1 1 concerns the parabolic problem: 



d t u(t,x) = d 2 xx {u m {t,x)), te]o,T], xeR, 
w(0,-) = S , 

where S is the Dirac measure at zero and u m denotes u|u| m_1 . It is well known that, 
for m > 1, there exists an exact solution to (1.2) , the so-called Barenblatt's density, see 
0. Its explicit formula is recalled for instance in 11341 Chapter 4] and more precisely in 
(4j Section 6.1]. Equation ( | 1 . 2 p is the classical porous medium equation. 

In this paper, we focus on ( |1.2| l when m e]0, 1[: the fast diffusion equation. In fact, 
an analogous Barenblatt type solution also exists in this case, see 11341 Chapter 4] and 
references therein; it is given by the expression 

U{t,x) =t- a (D + k\x\H- 2a ) T ^ T , (1.3) 

where 




m + 1 ' 2(m + l)m 



Equation ( |1.2| l is a particular case of the so-called generalized porous media type 
equation 

d t u(t,x) = dlJ{u{t,x)), te]0,T\, 
u(0,x) = uo(dx), x G R, 

where f3 : R — > R is a monotone non-decreasing function such that /3(0) =0 and uq is 
a finite measure. When f3{u) = u m , m e]0, 1[ and u = S , two difficulties arise: first, 
the coefficient f3 is of singular type since it is not locally Lipschitz, second, the initial 
condition is a measure. Another type of singular coefficient is /3(u) = H{u — u c )u, where 
H is a Heaviside function and u c > is some critical value, see e.g. [2]. Problem \1.2\ 
with m €]0, 1[ was studied by several authors. For a bounded integrable function as 
initial condition, the equation in \l.2\ is well-stated in the sense of distributions, as a by 
product of the classical papers 11101 on l|1.5|l with general monotonous coefficient (3. 
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When the initial data is locally integrable, existence was proved by H19B . 11 IP extended 
the validity of this result when u is a finite Radon measure in a bounded domain, 
02911 established existence when uq is a locally finite measure in the whole space. The 
Barenblatt's solution is an extended continuous solution as defined in H131 1140 : II 141 
Theorem 5.2] showed uniqueness in that class. If2~3 ' Theorem 3.6] showed existence in 
a bounded domain of solutions to the fast diffusion equation perturbed by a right-hand 
side source term, being a general finite and positive Borel measure. As far as we know, 
there is no uniqueness argument in the literature whenever the initial condition is a 
finite measure in the general sense of distributions. Among recent contributions, 01511 
investigated the large time behavior of solutions to l |1.2[ l. 

The present paper provides the probabilistic representation of the (Barenblatt's) 
solution of ( |1.2| l and exploits this fact in order to approach it via a Monte Carlo sim- 
ulation with an L 2 error around 10~ 3 . We make use of the probabilistic procedure 
developed in (U Section 4] and we compare it to the exact form of the solution U of 
fll.2| l which is given by the explicit formulae l [1.3| l-l fT~4"l ). The target of [4] was the case 
(3(u) = H(u — u c )u; in that paper those techniques were compared with a deterministic 
numerical analysis recently developed in II 1 2 11 which was very performing in that target 
case. At this stage, the implementation of the same deterministic method for the fast 
diffusion equation does not give satisfying results; this constitutes a further justification 
for the probabilistic representation. 

We define 

$(u) = \u\ m ^ 1 , ueE, m€]0,l[. 

The probabilistic representation of U consists in finding a suitable stochastic process Y 
such that the law of Y t hasU(t, ■) as density. Y will be a (weak) solution of the non-linear 
SDE 

f Y t = }V2$(U(s,Y s ))dW s , (16) 

{ U(t,-) = Law density of Y u Vfe]0,T], 

where W is a Brownian motion on some suitable filtered probability space (fi, J 7 , (J r t ) t >o, P)- 

To the best of our knowledge, the first author who considered a probabilistic rep- 
resentation of a solution of fll.5| ) was H. P. Jr. McKean ( 0260 ). particularly in relation 
with the so-called propagation of chaos. In his case /3 was smooth, but the equation 
also included a first order coefficient. From then on, literature steadily grew and nowa- 
days there is a vast amount of contributions to the subject, especially when the non- 
linearity is in the first order part, as e.g. in Burgers' equation. We refer the reader to 
the excellent survey papers 0330 and 0180 . A probabilistic interpretation of < |1 . 5| when 
j3(u) = to > 1, was provided for instance in [51- Recent developments related 

to chaos propagation when (3(u) — u 2 and f3{u) — u m ,m > 1 were proposed in 02811 
and Q17B . The probabilistic representation in the case of possibly discontinuous /3 was 
treated in [71 when j3 is non-degenerate and in 01 when fj is degenerate; the latter case 
includes the case j3(u) = H(u — u c )u. 

As a preamble to the probabilistic representation we make a simple, yet crucial 
observation. Let W be a standard Brownian motion. 

Proposition 1.1. Let ft : R R such that f3(u) = $ 2 (m).w, $ : R ->• R + and u be a 
probability real measure. 

Let Y be a solution to the problem 

Y t = Y + S ^{u{s,Y s ))dW s , 

u(t,-) = Law density of Y t , V t e}0,T}, 
it(0, •) = uo(dx). 
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Then u : [0, T] x R ->• R is solution to ( [TT5| |. 

Proof of the above result is based on the following lemma. 

Lemma 1.2. Let a : [0, T] x R — ► R + be measurable. Let (Y t ) be a process which solves 
the SDE 

Y t =Y 



■ J y/2a(s,Y a )dW s , te[0,T}. 



Consider the function t n- p(t, ■) from [0, T] to the space of finite real measures M(R), 
defined as pit, •) being the law ofY t . Then p is a solution, in the sense of distributions 
(see $Z72\ ), of 

d t u = dl x {au), t G]0,T], g 
u(0, •) = LawofFo- 

Proof of Lemma \l .21 This is a classical result, see for instance H32I Chapter 4]. The 
proof is based on an application of Ito's formula to <p{Y t ), ip e <S(R). □ 



Proof of Proposition We set a(s,y) — <& 2 (w(s, y)). We apply Lemma 1.2 setting 



p(t, dy) = u(t, y)dy, t e]0, T], and p(0, •) = u . □ 
When uq is the Dirac measure at zero and /3(u) — u m , with m €]§,![, Theorem 



5.7| states the converse of Proposition |1.1| providing a process Y which is the unique 
(weak) solution of \1.6) . The first step consists in reducing the proof of that Theorem 
to the proof of Proposition |5 . 3| where the Dirac measure, as initial condition of ( |1.2| l, is 
replaced by the function U(k, ■), < k < T. This corresponds to the shifted Barenblatt's 



solution along a time k, which will be denoted by U. Also, in this case Proposition 5.3 
provides an unique strong solution of the corresponding non-linear SDE. That reduction 
is possible through a weak convergence argument of the solutions given by Proposition 
|5.3| when k — > 0. The idea of the proof of Proposition |5.3| is the following. Let W be a 
standard Brownian motion and Y be a r.v. distributed as U(k, •); since is Lipschitz, 
the SDE 



Y t =Y a + f ${U( Sl Y s ))dW s , t G]0,T] 
Jo 



admits a unique strong solution. The marginal laws of (It) and U can be shown to be 
both solutions to ( | 1 . 6 P for a(s, y) = (U(s, y))" 1 ^ 1 ; that a will be denoted in the sequel by 
a. The leading argument of the proof is carried by Theorem |3. 1 | which states uniqueness 
for measure valued solutions of the Fokker-Planck type PDE ( |1.8[ ) under some Hypoth- 
esis(B). More precisely, to conclude that the marginal laws of (Y t ) and U coincide via 
Theorem|3.1| we show that they both verify the so-called Hypothesis(B2). In order to 



prove that for U, we will make use of Lemma 4.2 The verification of Hypothesis(B2) 
for the marginal laws of Y is more involved. It makes use of a small time (uniformly 
with respect to the initial condition) upper bound for the density of an inhomogeneous 
diffusion flow with linear growth (unbounded) smooth coefficients, even though the dif- 
fusion term is non-degenerate and all the derivatives are bounded. This is the object of 
Proposition [57T] the proof of which is based on an application of Malliavin calculus. In 
our opinion this result alone is of interest as we were not able to find it in the literature. 
When the paper was practically finished we discovered an interesting recent result of 
M. Pierre, presented in H201 Chapter 6], obtained independently. This result holds in 
dimension 1 when the coefficients are locally bounded, non-degenerate and the initial 
condition has a first moment. In this case, the hypothesis of type (B) is not needed. 
In particular it allows one to establish Proposition |5.3| but not Theorem |5.7| where the 
coefficients are not locally bounded on [0, T] x R. 
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The paper is organized as follows. Section [2] is devoted to basic notations. Section 
[3] is concentrated on Theorem |3.1| which concerns uniqueness for the deterministic, 
time inhomogeneous, Fokker-Planck type equation. Section [4] presents some proper- 
ties of the Barenblatt's solution U to < |1 -2| - The probabilistic representation of U is 
treated in Section [5] Proposition |5.1| performs small time density estimates for time- 
inhomogeneous diffusions, the proof of which is located in the Appendix. Finally, Sec- 
tion [6] is devoted to numerical experiments. 

2 Preliminaries 

We start with some basic analytical framework. In the whole paper d will be a strictly 
positive integer. If / : R d — >• R is a bounded function we will denote ||/||oo = sup |/(a;)|. 

x£R d 

By S(R d ) we denote the space of rapidly decreasing infinitely differentiable functions 
(p : R d — s- R, by S"(R d ) its dual (the space of tempered distributions). We denote by 
M(R d ) the set of finite Borel measures on R d . If x G R d , \x\ will denote the usual 
Euclidean norm. 

For e > 0, let K e be the Green's function of e — A, that is the kernel of the operator 

(e - A)" 1 : L 2 (R d ) -> H 2 (R d ) c L 2 (R d ). In particular, for all ip € L 2 (R d ), we have 

B e <p := (e - A)-VW = / K e (x - y) tp(y)dy. (2.1) 

For more information about the corresponding analysis, the reader can consult H31B . If 
tp € C 2 (R d ) p| S"(R d ), then (e— A)(p coincides with the classical associated PDE operator 
evaluated at ip. 

Definition 2.1. We will say that a function ip : [0,T]xl->l is non-degenerate if there 
is a constant c > such that tp > c„. 

Definition 2.2. We will say that a function ^:[0,T]xt->E has linear growth (with 
respect to the second variable) if there is a constant C such that \ip(-,x)\ < C(l + \x\), 
xeR. 

Definition 2.3. Let a : [0,T] x R d -> R+ be a Borel function, z° e M(R d ). A (weakly 
measurable) function z : [0, T] — > M(R d ) is said to be a solution in the sense of distribu- 
tions of 

dtz = A(az) 

with initial condition z(0, •) = z° if, for every t e [0, T] and 4> € <S(R), we have 

/ <fi(x)z(t, dx) = I cf>{x)z [dx) + ds Acj)(x)a(s,x)z(s,dx). (2.2) 
jR d Jn d Jo Jn d 

3 Uniqueness for the Fokker-Planck equation 

We now state the main result of the paper which concerns uniqueness for the Fokker- 
Planck type equation with measurable, time-dependent, (possibly degenerated and un- 
bounded) coefficients. It generalizes Theorem 3.8] where the coefficients were 
bounded and one-dimensional. 

The theorem below holds with two classes of hypotheses: (Bl), operating in the 
multidimensional case, and (B2), more specifically in the one-dimensional case. 

Theorem 3.1. Let a be a Borel nonnegative function on [0, T] x R d . Let z t : [0,T] — > 

M(R d ), i = 1,2, be continuous with respect to the weak topology on finite measures on 
M(R d ). Let z° be an element of M(R d ). Suppose that both Z\ and z 2 solve the problem 
d t z = A(az) in the sense of distributions with initial condition z(0, ■) = z°. 

Then z := (zi — z 2 )(t, •) is identically zero for every t under the following requirement. 
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Hypothesis (B). There is z e L} oc ([Q,T] x B d ) such that z(t, ■) admits z(t, •) as density 
for almost all t e [0, T]; z will still be denoted by z. Moreover, either (Bl) or (B2) below 
is fulfilled. 

(Bl) (i) J \z(t,x)\ 2 dtdx < +00, (ii) J \az\ 2 (t, x)dtdx < +00. 

[0,T]xK d [0,T]xK d 

(B2) We suppose d = 1. For every t > 0, we have 

(i) J \z(t, x)\ 2 dt dx < +00, (ii) J \az\(t, x) dtdx < +00, (hi) J \az\ 2 (t,x) dtdx < +00. 

[t ,T]xE [0,T]xE [t ,T]xR 

Remark 3.2. The weak continuity of z(t, •) and [7, Remark 3.10] imply that sup \\z(t,-)\\ var < 

te[o,T] 

+00, where || • \\ var denotes the total variation. In particular sup J Rd \z(t,x)\dx < +00. 

0<t<T 

Remark 3.3. 1. If a is bounded then the first item of Hypothesis(Bl) implies the 
second one. 

2. If a is non-degenerated, assumption (ii) of Hypothesis(Bl) implies assumption (i). 
Remark 3.4. Let d = l. 

1. If a is non-degenerate, the third assumption of Hypothesis(B2) implies the first 
one. 

2. If z(t, x) € L°°([t , T] x R) then the first item of Hypothesis(B2) is always verified. 

3. If a is bounded then assumption (ii) of Hypothesis(B2) is always verified by Remark 
|3.2| the first item of Hypothesis(B2) implies the third one. So Theorem |3.1| is a 
strict generalization of Theorem 3.8]. 

4. Let (z(t, 6 [0, T]) be the marginal law densities of a stochastic process Y solv- 
ing 

Y t =Y + f y/2a(s,Y s )dW S) 



with Y distributed as z° such that / |a;| 2 z (cfe) < +00. 

E. 

If y/a has linear growth, it is well known that supE(|Yi| 2 ) < +00; so 



/ Hs,x)z(s,x)\, 



J ds dx = E 

[0,T]xR 

Therefore assumption (ii) in Hypothesis(B2) is always fulfilled. 



t<T 
T 



a(s,Y s ) ds 



< +00. 



Proof of Theorem \3. 1 1 Let z\, z 2 be two solutions of l |2.2| l; we set z := z\ — z 2 . We 
evaluate, for every t e [0, T], the guantity 

ge(t) = \\z(t,-)\\\ e , 

where ||/||_ liE = ||(e - A)-*/IU». 

Similarly to the first part of the proof of [7, Theorem 3.8], assuming we can show 
that 

lims e (t) = 0, Vte[0,T], (3.1) 
we are able to prove that z(t) = for all t e]0, T}. We explain this fact. 
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Let t e]0, T]. We recall the notation BJ = (e - A)" 1 /, if / G i 2 (E. d ). Since z(t, •) € 
L 2 (R d ) then B e *(i, •) e iJ 2 (R d ) and so \7B £ z(t, •) e i? 1 (K, d ) d c L 2 (R d ) d . This gives 

<7 e (i) = / B e z(t,x)z(t,x)dx — e I (B £ z(t,x)) 2 dx — / B £ z(t,x)AB £ z(t,x)dx 
Jw JR d Jm d 



e j (B £ z(t,x)) 2 dx + / \VB £ z(t,x)\ 2 dx. 
JTs, d Jn d 



Since the two terms of the above sum are non-negative, if fl3.1| l holds, then ^/eB £ z(t, ■) — > 
(resp. \VB e z(t,-)\ -> 0) in L 2 (H d ) (resp. in L 2 (R d ) d ). So, for all t e]0,T], •) = 
eB e z(t, •) — AB £ z(t, •) — > 0, in the sense of distributions, as e goes to zero. Therefore 

Z EE 0. 

We proceed now with the proof of < |3. 1 1 . We have the following identities in the sense 
of distributions: 



z(t,-)= / A(az)(s,-)ds = / (A- e)(az)(s,-)ds + e (az)(s,-)ds, (3.2) 
Jo Jo Jo 

which implies 

B £ z{t,-) = - f (az)(s,-)ds + e [ B £ (az){s, -)ds. (3.3) 
Jo Jo 

Let 8 > and (<j>s) a sequence of mollifiers converging to the Dirac delta function at 
zero. We set z§(t, x) — f„ d z(t, y)4>s(x — y)dy, observing that zg € (L 1 f] L°°)([Q, T] x R d ). 
Moreover, ( |3.2| l gives 

zg(t, •) = / A(az) s (s, -)ds. 



We suppose now Hypothesis(Bl) (resp. (B2)). Let t„ = (resp. t > 0). By assumption 
(Bl)(ii) (resp. (B2)(iii)), we have A(az)s G L 2 ([i ,T] x Thus, z$ can be seen as a 
function belonging to C([t , T];L 2 (R d )). Besides, identities ( |3.2[ ) and < |3 . 3P lead to 



z s (t,-) = z s (t ,-)+ / (A-e)(az) 5 (s,-)ds + e (az) s (s, -)ds, (3.4) 

Jt Jt a 

B £ z s (t,-) = B £ z s (t ,-)- / (az) s (s,-)ds + e B £ (az) s {s, -)ds. (3.5) 

J t Jt Q 

Proceeding through integration by parts with values in L 2 (R d ), we get 

\\zs(t, W-l.e - IM*0, i, e = - 2 / ds < ^( S > ')> (<«)«(*, >L 2 



to 



(3.6) 



+ 2s / ds < (az)i(s, ■), B e z s (s, •) > L 2 



Then, letting 8 go to zero, using assumptions (Bl)(i)-(ii) (resp. (B2)(i) and (B2)(iii)) and 
Cauchy-Schwarz inequality, we obtain 

\\z(t, Oil — i ^ — lk(*o, •)ll-i, e = -2 / ds / a( S ,x)|z| 2 ( S ,a;)da; 

(3.7) 

+ 2e f ds < (az)(s,-),B £ z(s,-) > L i . 

J t 
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we have z(t , •) = 0. Using the inequality c t c 2 < Cl ~t° 2 , c lt c 2 € R and Cauchy-Schwarz, 



At this stage of the proof, we assume that Hypothesis(Bl) is satisfied. Since t = 0, 

2 +c 2 
2 

3.7| ) implies 

IK*,-)||-i, E < -2 f ds f (a\z\ 2 )(s,x)dx + e [ ds\\az(s, -)\\ 2 L2 + e [ ds\\B s z(s,-)\\ 2 L2 

J0 JR d JO Jo 

< ef ds\\az{s,-)f L 2+ [ ds\\z(s,-)f_ ltS) (3.8) 
Jo Jo 



because for / = z(s, •), we have 

We observe that the first integral of the right-hand side of ( |3.8[ l is finite by assumption 
(Bl)(ii). Gronwall's lemma, applied to H3.8|l, gives 



T 

T 



\z(t,-)\\^ l£ <ee / ds\\az(s, 



\l 2 - 



Letting e — >• 0, it follows that \\z{t, •)|| 2 _ 1 e = 0, Vt G [0, T]. This concludes the first part of 
the proof. 

We now suppose that Hypothesis (B 2) is satisfied, in particular d = l. By [7, Lemma 
2.2] we have 

sup2e\B e z(s,x)\ < y/e\\z(s,-)\\ var . 



Consequently l |3.7| l gives 

INC** Oll-i e — IN(*a, -)ll-i s < \/esup ||2;(*, -)IU«r / \az\(s,x)dsdx. (3.9) 

t<T J 

[t ,T]xR 

Besides, arguing like in the proof of [7, Theorem 3.8], we obtain that 

lira J|*(t o ,-)ll-i, s = 0. 
We first let t„ — >• in ( |3.9[ l, which implies 

\\z(t,-)f_ u < VEsap\\z(t,-)\\ var ( \az\{s,x)ds; (3.10) 



t<T 

[0,T]xR 

we remark that the right-hand side of ( |3 . 1 Op is finite by assumption (B2)(ii). Letting e 
go to zero, the proof of ( |3.1| l is finally established. 

□ 

4 Basic facts on the fast diffusion equation 

We go on providing some properties of the Barenblatt's solution U to \l.2\ when 
m e]0, 1[ and given by jl.3H1.4| ). 
Proposition 4.1. 

(i) U is a solution in the sense of distributions to l |1.2| l. In particular, for every <p G 
C§°(E), we have 

t 

' <p(x)U(t,x)dx = <p(0) + J ds j ' U m {s,x)ip"{x)dx. (4.1) 

R 
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(ii) J U(t, x)dx = 1, Vf > 0. In particular, for any t > 0, U(t, ■ ) is a probability density. 

R 

(Hi) The Dirac measure <5 is the initial trace ofU, in the sense that 



J i{x)U{t,x)dx ->7(0), as i^O, (4.2) 

B 



for every : R — > R, continuous and bounded. 



Proof of Proposition \4.1\ (i) This is a well known fact which can be established by in- 
spection. 

(ii) For M > 1, we consider a sequence of smooth functions (f M ), such that 

f = 0, if \x\ > M + 1; 

f M (x)l < 1, if \x\ £ [M, M + 1]; 

i = 1, if \x\ < M. 

By l |4.1[ l we have 



tp M {x)U{t,x)dx = 1 + / ds / U m (s,x)((p M ) "(x)dx. (4.3) 

E OR 

Letting M — > +oo, by Lebesgue's dominated convergence theorem, the left-hand side of 
( |4.3[ l converges to J R U(t,x)dx. The integral on the right-hand side of ( |4.3| l is bounded 
by 

t Af+l t T 

C J ds J U m (s,x)dx < C J s- am (D + feM 2 s" 2a ) ^ ds < m / s^ds. 

OA/ ^ 

The last integral on the right is finite as ^ > 0, for every m €]0, 1[. Therefore the 
integral in the right-hand side of ( |4.3[ l goes to zero as M — > +oo. This concludes the 
proof of the second item of Proposition |4.1| 

(hi) (|4.2[l follows by elementary changes of variables. □ 



Note that the second item of Proposition |4.1| determines the explicit expression of 
the constant D. 

Lemma 4.2. 

(i) Suppose that ^ < m < 1 . Then there is p > 2 and a constant C p (depending onT) 
such that for < s < I <T 

J dtdx(U(t,x)) Si ^ +1 <C p {i- s). (4.4) 

]s,!]xE 



(ii) In particular, taking p = 2 in H4.4| l, we get 

dtdx{U{t,x)) m < +oo, (4.5) 



]0,T]xR 

again when m belongs to ] i , 1 [. 
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(in) If j < m < 1 , 



dtdx(U(t,x)) 2m < +00. (4.6) 



]0,T]xR 

(iv) If m belongs to } 3 . 1 [ , then 



Vk>0, J \x\ l U{n,x)dx <+oo. (4.7) 



Proof of Lemma \4^2l 

(i) Using l|1.3|l, we have 



]a,«]xE ]s,(]xR 

Then, setting y = t~ a xy j^, we get 



p 1 

dtdx. 



{U(t,x)) lA ^ +1 dtdx^ D 2 jt^^dt [(l + y 2 )%-^dy 

]s,f]xE s 

< 2 — Jl r^-^^-a) / {l + y 2 )i-^dy. 

Vk 

The last integral is finite if (p+ 1)(1 — m) < 2. This implies l |4.4| ). 

(ii) is a particular case of (i) and (iii) follows by similar arguments as for the proof of 

(i). 

(iv) Now we assume that m s] |, 1[. For n > we have 

J \x\ i U(K,x)dx = I> p /2 °'' k 4 " J \y\*(l + y 2 )-^dy, (4.8) 

R R 

where this last equality was obtained setting y = n~ a x\J^. Clearly, since m e]|, 1[, 
the integral in the right-hand side of (|4.8[l is finite. Therefore (|4.7|l is fulfilled. □ 



Let /c e]0, T}. Given u : [0, T] x R -> R we associate 

u(t,x)=u(t + K,x), (t, x) G [0, T - k] x R. (4.9) 

In particular we have 

U(t,x)=U(t + K,x). (4.10) 
Moreover, for every 2; e R, we denote 

u 0:H (x) =U(k,x). (4.11) 

Remark 4.3. Function U solves the problem 



d t u = d 2 xx {u m ), 
u(Q, •) = w ,k • 



(4.12) 
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5 The probabilistic representation of the fast diffusion equation 

We are now interested in a non-linear stochastic differential equation rendering the 
probabilistic representation related to < |1 -2| and given by l |1.6[ l. Suppose for a moment 
that Y is a random variable distributed according to S , so Y„ = a.s. We recall that, 
if there exists a process Y being a solution in law of ( |1.6| l, then Proposition 1 1 . 1 1 implies 
that u solves \1.2) in the sense of distributions. 

In this subsection we shall prove existence and uniqueness of solutions in law for 
HI. 6} , In this respect we first state a tool, given by Proposition |5 . 1 1 below, concerning 
the existence of an upper bound for the marginal law densities of the solution Y of 
an inhomogeneous SDE with unbounded coefficients. This result has an independent 
interest. 

Proposition 5.1. Let <r, b : [0, T] x 1R — >• 1R be continuous (not necessarily bounded) 
functions such that cr(t, ■), b(t, •) are smooth with bounded derivatives of orders greater 
or equal than one. a is supposed to be non-degenerate. 
Letx G R and Y t = (yT°)te[o,T] be the solution of 

Y t = x + [ a(r,Y r )dW r + f b(r,Y r )dr. (5.1) 
Jo Jo 

Then, for every s > 0, the law ofY s admits a density denoted p s (x ,-). 
Moreover, we have 

Ps (x , x) < -p (1 + K| 4 ) , V(s, x) G]0, T] x R, (5.2) 

where K is a constant which depends on \\a'\\ x , ||6'||oo and T but not on x . 
Remark 5.2. 1. The proof of Proposition [5~T] above is given in Appendix |7.l| 

2. If a and b is bounded, the classical Aronson's estimates implies that ( |5.2[ l holds 
even without the |a; Q | 4 multiplicative term. If a and b are unbounded, 1 1 1 provides 
an adaptation of Aronson's estimates; unfortunately they first considered time- 
homogeneous coefficients, and also their result does not imply l |5.2[ l. 

3. If a and b have polynomial growth and are time-homogeneous, various estimates 
are given in 02511 . However the behavior is of type 0(t~i) instead of 0{t~*) when 

Let Y K be a random variable distributed according to u 0:K . We are interested in the 
following result. 

Proposition 5.3. Assume that m 1[. Let B be a classical Brownian motion inde- 
pendent ofY K . Then there exists a unique (strong) solution Y = {Y t ) te [ 0jT _ K ] of 



Y t = Y K + j $(U( S ,Y s ))dB s , 

_ o 

U(t,-) = Law density of Y t , Vte[0,T 

W(0,-) = u . K . 



(5.3) 



In particular pathwise uniqueness holds. 

Corollary 5.4. Let W be a classical Brownian motion independent ofY K . Therefore 
there is a unique (strong) solution Y K = {Y t K ) telK/r] of 

Y t K = Y K + mU(s,Y«))dW s , 
U(t,-) = Law density of Y t K , V t e [k,T], (5 ' 4) 

U(K, ■) = !t„, K . 
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Proof of CoroUary \5.4\ We start with the proof of uniqueness. Let n > 0. We consider 
two solutions Y K ' 1 and Y K ' 2 of $5A) , we set Yf = Y t K / K , Wt G [0,T - k], i = 1,2 and 
B t = Wt +fc - Wt, Vf G [0,T - «]. Clearly and if solve (53) . Therefore, using 
Proposition |5.3| we deduce uniqueness for problem ( |5.4| l. Existence follows by similar 
arguments. □ 



Proof of Proposition \5.3\ Let W be a classical Brownian motion on some filtered prob- 
ability space. Given the function U, defined in ( |4. 1 0| >, we construct below a unique 
process Y strong solution of 

t 

Y t = Y + J <P(U(s,Ys))dW s . (5.5) 
o 

From ( |4. 1 Q[ >, for every (s, y) e [0, T — k] x EL, we have 

$(fii(s,y)) = y/2a(s,y), 

where 

a(s,y) = (s + n) ail - m) (D + k\y\ 2 (s + n)- 2a ). (5.6) 

In fact, is continuous, smooth with respect to the space parameter and all the 

space derivatives of order greater or equal than one are bounded; in particular <&(U) is 
Lipschitz and it has linear growth. Therefore fl5.5| ) admits a strong solution. 



By Lemma 1.2 the function t H> p(t, •) from [0,T— k] to VW(EL), where p(t, •) is the law 
of Y tl is a solution to 

dtp = dl x (ap), 
P(0, ■) = u 0tK . 

To conclude it remains to prove that U(t, y)dy is the law of Y t , Vt G [0, T —k]; in particular 
the law of the r.v. Y t admits a density. For this we will apply Theorem |3.l| for which we 
need to check the validity of Hypothesis(B2) when a = a and for z := z\ — z 2 , where 
z\\= p and z 2 -=U. By additivity this will be of course fulfilled if we prove it separately 
for z := p and z :=U, which are both solutions to ( |5.7[ l. 

Since a is non-degenerate, by Remark |3.4{ l), we only need to check items (ii) and 
(hi) of the mentioned Hypothesis(B2). On one hand, since a(s,y) = U (s,y), z :=U 
verifies Hypothesis(B2) because of items (ii) and (hi) of Lemma |4~2j On the other hand, 



since \Ja has linear growth, by Remark 3.4 (4) p fulfills item (ii) of Hypothesis(B2). 
y Lemma |5.5| below 
implies that U = p. 



Moreove r, by Lemma [5.5| below, p also verifies item (hi) of Hypothesis(B2). Finally 
Theorem 

□ 

Lemma 5.5. Let : [0,T]xE-> R + , continuous (not necessarily bounded) such that 
ip(t, •) is smooth with bounded derivatives of orders greater or equal than one. We also 
suppose ij> to be non-degenerate. 

We consider a stochastic process X = {X t )te[o.T} strong solution of the SDE 

t 

X t =X + J iP(s,X s )dW s , (5.8) 
o 

where X is a random variable distributed according to u K defined in l |4.11[ l with m G 

]f.l[- 

For t g]0, T] the law of X t has a density v(t, ■) such that (tp 2 v)(t,x) belongs to 
L 2 ([t B ,T] x EL), for every t„ > 0. 
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Proof of Lemma [5751 

If X a — x , where x is a real number, then Proposition |5.1| implies that, for every 
t e]0, T], the law of X t admits a density p t (x ,-). Consequently, if the law of X is 
uo. K (x)dx, for every t g]0, T], the law of X t has a density given by 



v(t,x)= / u 0<K (x )p t (x ,x)dx . 
By (| 5 . 2 p in Proposition |5.1| it follows 



B 



K 

sup Pt(x ,x) < K (l + |a; | 4 ), where if = —=. 

(t,x)e[t ,T]xR V*o 



Using l|5.9[l we get 



ifi := sup \v(t,x)\<Ko I (1 + |x | 4 )W(k, x )dx < oo; 

(t,i)e[( ,T]xK 



(5.9) 



(5.10) 



the latter inequality is valid because of ( |4.7[ l in Lemma |4.2| In the sequel of the proof, 
the constants K 2 ,K 3 , K± will only depend on t , T and tj). Furthermore 



{{il?v){t,x)j 2 dtdx < sup \v(t,x)\E 

[t ,T]xK 

Since ip has linear growth, this expression is bounded by 



i>\t,X t )dt 









(i + /b 


sup |X t | 4 




\ 


te[o,r] 





(5.11) 



( |5.11[ l follows because of ( |5 . 1 OP . Besides, by Burkholder-Davis-Gundy and Jensen's in- 
equalities, taking into account the linear growth of ip, it follows that 









T 






E 


sup |X f | 4 


< K 3 [|X | 4 ] - 


f / E 


sup \X S \ 4 






te[o,T] 





se[o,T] 





Then, by Gronwall's lemma, there is another constant K 4 such that 



E 



sup \X t \ 
te[o,T] 



<K A \1 



\x \ 4 U{K,x )dx 



Finally ( |5.11| l, l |5.12| l and ( | 5 . 1 Op allow us to conclude the proof. 



(5.12) 



□ 



We are now ready to provide the probabilistic representation related to function U 
which in fact is only a solution in law of l |1.6| ). 

Definition 5.6. We say that ( ] 1 . 6 p admits a weak (in law) solution if there is a probability 
space (Q : J 7 , P), a Brownian motion (W t )t>o and a process (Y t ) t>Q such that the system 
d 1 - 6 p holds. ( |1.6| l admits uniqueness in law if, given (W 1 , Y 1 ), (W 2 , Y 2 ) solving d 1 . 6 P on 
some related probability space, it follows that Y 1 and Y 2 have the same law. 

Theorem 5.7. Assume that m 1[. Then there is a unique weak solution (in law) Y 



of problem d 1 . 6 



Elect ron. J. Probab. (2012), no. 0, 1- [28| 



ejp.ejpecp.org 



14 



Fokker-Planck and Fast Diffusion 



Remark 5.8. Indeed the assumption on m 1[ is only required for the application of 
The arguments following the present proof only use m > | . 



Theorem 



3.1 



Proof of Theorem \5. 71 First we start with the existence of a weak solution for ( |1 .6| . 
Let U be again the (Barenblatt's) solution of \1.2) . We consider the solution (Y t K ) te [ K) T) 
provided by Corollary 5.4 extended to [0, n], setting Y t K = Y K , t 6 [0, k]. We prove that 
the laws of processes Y K are tight. For this we implement the classical Kolmogorov's 
criterion, see [22, Section 2.4, Problem 4.11 ]. We will show the existence of p > 2 such 
that 

E[|r t K -y s K | p ] <C p \t-a\i, VMe[0,T], (5.13) 



where C p will stand for a constant (not always the same), depending on p and T but not 
on k. Let s, i e]0, T]. Let p > 2. By Burkholder-Davis-Gundy inequality we obtain 



K[\Y t " 



y s T] < c p e 



$ 2 (W(r,y r K ))dr 



Then, using Jensen's inequality and the fact that U(r, ■ 
we get 



is the law density of Y r K , r > k, 



We have 



E[\Y t K -Y s K \*] <C p \t-s\^ I dr J &>(U(r,y))U(r,y)dy. 



dr / &{U{r,y))U{r,y)dy= / dr / dy{U(r,y)) 



(5.14) 



and, by Lemma [472] (i), the result follows. 

Consequently there is a subsequence Y n := Y Kn converging in law (as C([0, T])— valued 
random elements) to some process Y. Let P n be the corresponding laws on the canon- 
ical space f2 = C([0,T]) equipped with the Borel cr-field. Y will denote the canonical 
process Y t (uj) — uj(t). Let P be the weak limit of (P n ). 

1 ) We first observe that the marginal laws of Y under P n converge to the marginal 
law of Y under P. Let t e]0, T]. If the sequence (k„) is lower than t, then the law of Y t 
under P n equals the constant law U(t, x)dx. Consequently, for every t e]0, T], the law of 
Yt under P is U(t, x)dx. 

2) We now prove that Y is a (weak) solution of ( |1.6| l, under P. By similar arguments 
as for the classical stochastic differential equations, see H321 Chapter 6], it is enough 
to prove that Y (under P) fulfills the martingale problem i.e., for every / e (R), the 
process 



(mp) /(y t )-/(o) 



f"{Y s )^ 2 {U{s,Y s ))ds, 



is an (Tg) -martingale, where (J" s ) is the canonical filtration associated with Y. Cf(R) 
stands for the set {/ e C 2 (R)|/, /', /" bounded }. Let E (resp. E") be the expectation 
operator with respect to P (resp. P n ). Let s, t € [0, T] with s < t and i? = R(Y r , r < s) be 
an T B — measurable, bounded and continuous (on C([0, T])) random variable. In order to 
show the martingale property (MP) of Y, we have to prove that 



E 



f(Yt) - f(Y s ) 



f"(Y r )$ 2 (U(r,Y r ))dr\R 



0, /eC 2 (R). 



(5.15) 
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We first consider the case when s > 0. There is n > n , such that n n < s. Let / e C 2 (R); 
since (Y s ) s > Kn , under P n , are still martingales we have 



E" 



f(Y f ) - f(Y s ) - i / f"{Y r )<S> 2 {U{r,Y r ))dr ) R 



= 0. 



(5.16) 



We are able to prove that l |5.15[ l follows from ( |5. 1 6| >. Let e > and N > such that 

t 



Jdr J U m (r,y)dy<e, 



(5.17) 



tol>#-i} 



where C is the linear growth constant of <f> 2 o U in the sense of Definition 2.2 In order 
to conclude, passing to the limit in l |5.16| l, we will only have to show that 



lim E" [F(Y)} - E [F(Y)} = 0, 

71— » + 00 



(5.18) 



where F(£) = J dr<S> 2 (U(r, i(r))f"(l(r))R(£{£),£, < s), F : C([0, T]) -> E being continuous 

s 

but not bounded. The left-hand side of l |5.18| l equals 

E™ [F(Y) - F N {Y)] + E" - E [^(Y)] + E [F N (Y) - F(Y)] 



(5.19) 



where 



£ 1 {n,N)+£ 2 {n,N)+£ 3 {n,N), 



F N (£) = / dr {<5> 2 {U(r, *(r)) A JV) /" (e(r))R(i(£), £ < «). 



Since k„ < s, for TV large enough, we get 



^(n^l^plUllHoo/dr / (<S> 2 (U(r,y)-N)U(r,y)dy 

s {$ 2 {U(r,y))>N} 
t 

< \\R\U\f Woo J dr J U m (r,y)dy<e\\R\\ 00 \\f"\\ 00 , 



(5.20) 



{|s/|>f-i} 



taking into account ( |5.17[ l and the second item of Lemma |4.2| For fixed N, chosen in 
( |5.17[ ), we have lim„^ +oc £2(1%, N) = 0, since F N is bounded and continuous. Again, 
since the law density under P of Y t , t > s , is U(t, •), similarly as for ( |5.20[ l, we obtain 
\£a(n,N)\ < £||i?||oo||/"||oo- Finally, coming back to l [5~T9l l, it follows 

lim sup |E" [F(Y)]-V[F(Y)}\ < 2e\\R\\ 00 \\f"\\ ooi 

n— >-+oo 



since e > is arbitrary, l |5.18[ > is established. So ( ]5.15| ) is verified for s > 0. 
3) Now, we consider the case when s = 0. We first prove that 



± 

E J <$> 2 {U{r,Y r ))dr 

.0 



< +OO. 



(5.21) 
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By item 1) of this proof, the law of Y r , r > admits U(r, •) as density. Consequently the 
left-hand side of fl5.21| l gives 



dr j &(U{r,y))U(r,y)dy = j dr j U m (r,y)dy, 

]0,T] R OR 



which is finite by the second item of Lemma |4~2j Coming back to l |5.15| l, we can now let s 
go to zero. Since Y is continuous and / is bounded, we clearly have lim^o E [f(Y s )R] = 
E[f(Y )R]. Moreover 



limE 

s->0 



f"(Y r )<S> 2 (U(r,Y r ))dr \R 











)« 


= E 




;/ 








\0 



f(Y r )<S> 2 (U(r,Y r ))dr \R 



using Lebesgue's dominated convergence theorem and l |5.21[ l. Consequently we obtain 



E 



f(Y t ) - f(Y ) - - I f"(Y r )^(U(r,Y r ))dr | R 



= 0. 



(5.22) 



It remains to show that Y = a.s. This follows because Y t — > Y a.s., and also in law (to 
S ), by the third item of Proposition |4.1| Finally we have shown that the limiting process 
Y verifies (MP), which proves the existence of solutions to | |1.6| ). 

4) We now prove uniqueness. Since U is fixed, only uniqueness for the first line of 
equation (1.6| > has to be established. Let (5?)te[o,Ti- i = 1,2, be two solutions. In order 
to show that the laws of Y 1 and Y 2 are identical, according to 0211 Lemma 2.5], we 
will verify that their finite marginal distributions are the same. For this we consider 
= t < t x < ... < t N = T. Let < k < i x . Obviously we have Y t l a = a.s., in 
the corresponding probability space, Vi € {1,2}. Both restrictions F^ut] and F 2 |[ K! t] 
verify ( |5.4[ l. Since that equation admits pathwise uniqueness, it also admits uniqueness 
in law by Yamada-Watanabe theorem. Consequently F^r^T] and F 2 |[ k ,t] have the same 
law and in conclusion the law of (Y^ , . . . , Y^ ) coincides with the law of ( Y 2 x , . . . , Y 2 N ), 



thus, the law of (Y^ , 



Y 1 
• 1 t N 



coincides with the law of (Y 2 g , . . . , Y 2 ^ 



□ 



6 Numerical experiments 

In order to avoid singularity problems due to the initial condition being a Dirac delta 
function, we will consider a time translation of U, denoted v, and defined by 

«(*,•) =U(t + l,-), Vie [0,T]. 

v still solves equation ( |1.2| l, for m €]0, 1[, but with now a smooth initial data given by 

v (x) =U(l,x), Vxel. (6.1) 

Indeed, we have the following formula 

_ i 

v{t,x) = {t + iy a [D+~k\x\ 2 {t + l)- 2a ^j 1_m , (6.2) 

where a, k and D are still given by \1A) . 

We now wish to compare the exact solution of problem (L2) to a numerical proba- 
bilistic solution. In fact, in order to perform such approximated solutions, we use the 
algorithm described in Sections 4 of [4] (implemented in Matlab). We focus on the case 
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171=1- 

Simulation experiments: we compute the numerical solution over the time-space 
grid [0, 1.5] x [—15, 15]. We use n — 50000 particles and a time step At = 2x 10~ 4 . Figures 
[T](a)-(b)-(c)-(d) / display the exact and the numerical solutions at times t = 0, t = 0.5, 
t = 1 and t = T = 1.5, respectively. The exact solution for the fast diffusion eguation 
fll.2| l, given in H6.2D , is depicted by solid lines. Besides, Figure [T](e) describes the time 
evolution of the discrete L 2 error on the time interval [0, 1.5]. 




.2 
O.I 5 
0.1 

o.os 




10 15 



Figure 1: Numerical (dashed line) and exact solutions (solid line) values at t=0 (a), 
t=0.5 (b), t=l (c) and t=1.5 (d). The evolution of the L 2 error over the time interval 
[0,1.5] (e). 



7 Appendix 

7 . 1 Proof of Proposition |5.1| 

We start with some notations for the Malliavin calculus. The set D°° represents the 
classical Sobolev-Malliavin space of smooth test random variables. D 1 ' 2 is defined in the 
lines after 11271 Lemma 1.2.2] and L 1,2 is introduced in [27, Definition 1.3.2]. See also 
02411 for a complete monograph on Malliavin calculus. We state a preliminary result. 

Proposition 7.1. Let N be a non-negative random variable. Suppose, for every p > 1, 
the existence of constants C(p) and e (p) such that 

P(N<e)<C(p).eP+\ Vee]0,e o (p)]. (7.1) 

Then, for every p > 1, 

E(iV- p ) < e (p)-C(p + 1) + e (p)-PP(N > e (p)). (7.2) 

Proof of Proposition \71\ Let p > 1 and t {p) > 0. Setting F(x) = P(N < x), x e R+, we 
have 

E(JV-p) = h + h, (7.3) 
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where 

«o(p) +00 



h= J x~ p dF(x) and I 2 = J x- p dF(a 



eo(p) 

( |7.1[ l implies that 7i and i"2 are well-defined. Indeed, on one hand, applying integration 
by parts on I±, we get 

e o(p) 

h = [x- p F{x)Y° ip) + p J x- p - 1 F(x)dx; 



moreover, there is a constant C(p) such that 

h<(p+l)e (p)C(p). (7.4) 
On the other hand, again ( |7.1| l says that 

h<e (p)- p (l-F(e (p))). (7.5) 
Consequently, using (|7.4|l and (|7.5[l and coming back to (|7.3|l, l|7.2|l is established. □ 



Pi(x ,y) < -(l + |a; | 4 ), yyeR,\G}0,Vf], (7.6) 



Proof of Proposition \5.1\ In this proof a' (resp. b') stands for d x a (resp. d x b). Let Y = 
(Yt )te[Q.T\> be the solution of ( |5.1[ l. According to 0271 Theorem 2.2.2] we have Y s £ D°°, 
Vs £ [0,T]. Let s > 0. Since er is non-degenerate, by 0271 Theorem 2.3.1], the law of Y s 
admits a density that we denote by p s {x , •). 

The second step consists in a re-scaling, transforming the time s into a noise multi- 
plicative parameter A; we set A = yfs. Indeed, (Y t ) is distributed as (Ft ), where 



Y t x = x + A J a{r\ 2 ,Y r x )dW r + A 2 J b{r\ 2 .Y x )dr. 


In particular, Y s ~ Y x . By previous arguments, for every t > 0, Y t x e D°° and its law 
admits a density denoted by p x (x , ■). Our aim consists in showing the existence of a 
constant K such that 

K 
~\ 

where K is a constant which does not depend on x and A. In fact, we will prove that, 
for every A e]0, VT], 

sup p x t {x ,y) < ^(l + \x \ 4 ). (7.7) 
yeK,te]o,i] A 

We set Zl = Yt ^ Xo , t e [0, 1], so that the density of Z t A fulfills q£(z) = Xp x (x , Xz + 
x ), (t, z) £ [0, 1] x R. In fact, we will have attained ( |7.7| l, if we show 

sup q X (z)<K{l + \x f), te}0 7 l}- (7.8) 

zer.,\e]o,Vt) 

We express the equation fulfilled by Z; it yields 

t t 
Z x = J a x (r,Z x )dW r + J b x {r,Z x )dr, (7.9) 

where, for every (r, z) £ [0, 1] x R, we set 

cr x (r, z) — a(r\ 2 , Xz + x ), and b x (r, z) = Xb(rX 2 ,Xz + x ). 

At this stage we state the following lemma. 
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Lemma 7.2. For A e]0, 1], we shorten by Z := (2 ( A ) fg [ 0J ], the solution of ( |7.9| l. For 

every j > 1, we have 



sup E 



sup \Z t \ 
te[o,i] 



<C(l + |a: o n 



(7.10) 



where C is a constant depending on ||c/||oo, ||&'||oo and ?\ but not on x . 
Remark 7.3. 



1. For simplicity, in the whole proof of Proposition JETT1 we will set T = 1. 

2. Since there is no more ambiguity, we will use again the letter s in the considered 
integrals. 

Proof of Lemma [7.2\ Let A €]0, 1] and 7 > 1. In the proof C\ is a constant depending 
on T, and C 2 , C3 depend on T, ||ct'||oo and ||6'||oo- Using Burkholder-Davis-Gundy and 
Jensen's inegualities, we get 

<Ci Qf E[|ct(sA 2 ,AZ s +x )p}ds + J^ I][\\b(s\ 2 ,XZ s +x )p}ds 

Since <r' and 6' are bounded, a and have linear growth. Therefore, previous expression 
is bounded by 



E 


sup \z p y 


Ml 




pe[o,t] 





C 2 (l + A 7 ) 1 + |x 7 + A 7 / E 









/ E 


sup |Z p | 7 




/o 


_ P e[o,s] 





Since A e]0, 1] we obtain 
E 



sup |Z p | 7 
pe[o,t] 



Consequently, using Gronwall's lemma, the result follows. 



<C 3 [l + \x P + / E 









/ E 


sup |Z p | 7 


dsj 


Iq 


_ P e[o,s] 





□ 



Now, in order to perform (|7.8[l, we make use of Malliavin calculus for deriving ex- 



t t 



pression l |7.9[ l. Omitting A in the notation Zf , we get 

t t 
D r Z t = ct(VA 2 , AZ r +a; )l [ra] (t)+A f a' (s\ 2 , XZ s +x )D r Z s dW s +X 2 [ V (sX 2 , XZ s +x )D r Z s ds. 

Consequently 

D r Z t = a{rX 2 , XZ r + x )£ I A J a 1 {sX 2 , XZ S + x )dW s + X 2 J b' (sA 2 , XZ S + x )ds I , r /. 

\ r r / 

where £ (S) denotes the Doleans exponential of the continuous semi-martingale 

t t 
St=xJ a'(sX 2 ,XZ s +x )dW s + X 2 J b'(sX 2 ,XZ s +x )ds, te[0,l}. 

r r 

We recall that, £ t {S) = exp(£ t — §[5]t). Consequently, for fixed t e]0, 1], we have 

t 

< DZ t ,DZ t >= J a 2 {rX 2 ,XZ r + x )£ 2 (X;r,t)dr, (7.11) 
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where 



£(A; r, t)=explxj a'(sX 2 ,XZ s + x Q )dW s + A 2 /(&' - ^)(sX 2 , XZ S + x )ds J . (7.1 



2) 



We set, for every s < i, G(A, s) = — r)z s DZ > ' In view °^ tne application of Il27l 
Proposition 2.1.1], which implies the useful expression ( |7.25| l for the density of Z t , 
we will need to show that G(A, •) belongs to the domain of the divergence operator S, 
denoted by Dora S. It will be the case if G(X, •) £ l}' 2 (H) with H = L 2 ([0,T}). In fact, 
by the lines after |27, Definition 1.3.2], we know that E 1 ' 2 c Dom 5. Since Z t £ D°°, we 
can deduce that < jjz ^ DZ > Be l° n S s t0 provided that we prove 



eL p (n), Vp>i. 

<DZ u DZ t > V h 1 - 1 



(7.13) 



see E3 Lemma 2.1.6]. Since D°° is an algebra, G(A,s) £ D™, for s e]0,T], and so 
G(X,s) £ D 12 . < [7TT3] > will be the object of Proposition [73] According to (271 Defini- 
tion 1.3.2], to affirm that G(A, •) belongs to E 1 ' 2 it remains to show the existence of a 
measurable version of D Sl G(X, s), (si, s) <G [0, i] 2 , such that 



E 



(D Sl G(X,s)) 2 d Sl ds 



[0,t]- 



< +oo. 



(7.14) 



We first state the following Lemma. 

Lemma 7.4. For every q > 1, there exists a constant C (q) such that 

sup E[(£(X;r,t)y] <C (q). 

0<r<t<l,AG]0,l] 

Proof of Lemma VTM Let A e]0, 1] and g > 1. For fixed < r < t < 1, ([7TT2]) gives 



(7.15) 



£ 9 (A;r,i) = M(A;r,i,g)exp ( A (<j g) / (a') 2 (pA 2 , AZ P + x )dp + q\ 2 I b'( P X 2 ,XZ p + x )dp 



where 



M(X;r,t,q) = exp I / Agcr'OA 2 , AZ P + a^dWp - - J {qXa') 2 (pX 2 , XZ P + x )dp 



(7.16) 



In fact, since a' is bounded, the stochastic exponential M(X;r,t,q) verifies Novikov's 
condition; therefore it is a martingale. So E(M(A; r, t, q)) = 1. In addition, since b' is also 
bounded and A e]0, 1], wegetE[(£ (X;r,t))i] < C (q), where C (q) = exp (2{q 2 - q)\\a'\\1 + 2q\\b'\\ 
Conseguently l |7.15| l is established. □ 

Proposition 7.5. There is a constant C (not depending on x ) such that 



sup E[(< DZ t ,DZ t >y p ] <C, Vp > 1. 
(t,A)e]o,rj 2 

Proof of Proposition \7.5\ Let t e]0, 1] fixed, e = where c„ is a non-degeneracy con- 



stant of o~ in the sense of Definition 2.1 Consider e €]0, e Q [, we set N :— N x =< 
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DZ£,DZ* >, where we recall that < DZ t , DZ t > appears in ( |7.11[ l and l |7.12| l. Accord- 
ing to Proposition |7.1| we have to evaluate P(N < e). Taking into account the lower 
bound of a we have 



P(N<e) < p(J dr£ 2 (X;r,t) < ^\ 



(7.17) 



< P 



< P 



dr£ 2 (\-r,t)\ < 



dr 



£ 2 {X;r,t)dr\ > a — 



By the inverse triangle ineguality of the L'Qt — — ,t])-nonn we get 



P(N<e)<P[f (1 -£(X;r,t)) 2 dr > — | 
\Jt-M ~ c ) 



Let p > 1. By Chebyshev's ineguality this is lower than 



Co\ p+1 



(?) 



E 



(1 -£(X;r,t)) 2 dr 



Then, using Jensen's ineguality we get 

-1 r t 



P{N < e) < 4 p 



E[(l-£(X;r,t)) 2{p+1) ]dr. 



(7.18) 



Furthermore l |7.12[ l implies that £ (A; r, t) solves 

£(A; r,t) = 1 + X J £(X; r, s)a\sX 2 ,XZ s + x )dW s + A 2 / t\ A: r. * )//( *X\ AZ, -h .«•„ )<l*. 
Thus 



E [(£(A;r,t) - l) 2(p+1 >] < 2 2(p+1) E 



A / £{X;r,s)cr'(sX 2 , XZ S + x )dW s 



. 2 2 ("+ 1 >E 



A 2 / £(X;r,s)b'(sX 2 ,XZ s + x )ds 



2(p + l)- 



(7.19) 



On one hand, using Jensen's ineguality and A e]0, 1], we obtain 



E 



A 2 / £(X;r,s)b'(sX 2 ,XZ s + x )ds 



< \\b'\\Jt-r\ 2 ^ I E[(£(A;r, S )) 2(p+1 >]d S . 

(7.20) 



On the other hand, by Burkholder-Davis-Gundy ineguality we get 



E 



A / £{X;r lS y(sX 2 ,XZ s +x )dW s 



< E 



A 2 / £ 2 (X- 1 r,s)(a') 2 (sX 2 ,XZ s + x Q )ds 



p+r 
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Applying again Jensen's ineguality gives 



E 



< h'Ut-r\P I E[(£(\;r,s)Y^]d S . 

(7.21) 



A / £(X;r,s)a'{sX 2 ,XZ s + x )dW i 

r 

Therefore ( [7720) , %FT\) and (7719) lead to 

t 

E[(£(A;r,t)-l) 2<p+1) ] < C(T, ||a'|U ||&'|U)|i - rf / E [(£(A; r, s)) 2(p+1) ] ds. (7.22) 



By Lemma 7.4 and l |7.22[ l, there is a constant Co(2(p + 1)) such that 

E[(£(A;r,t)-l) 2(p+1) ] < d(T, ||a'|U, ||6'|U)C (2(p + 1)). (7.23) 
Then, coming back to l |7.18| l and using ( [7.23| l, we obtain 

Vee]0,e o ], P(N<e) <C{p)e» +1 , (7.24) 

where C(p) = ^""^^(p + ^KM? 1 , II^IU. II&'IL) . Finally using Proposition the 
result follows. □ 

We go on with the proof of Proposition |5.1| taking into account the considerations 

=]0,1], 

(7.25) 



before Lemma 7.4 In fact [27, Proposition 2.1.1] allows us to express, for fixed t €]0, 1], 

q?(z)=V[t {Zi>z} 5(G(X,-))] ; 



using Cauchy-Schwarz inequality it implies that 



q}{z)< x E |*(G(V))| 



(7.26) 



According to (1.48) in [27], l |7.26[ l implies 





" t 






q?(z) < ME 


J G 2 (X,s)ds 


+ E 


J (D Sl G(X,s)) 2 d Sl ds 




.0 







(7.27) 



Now we state a result that estimates the two terms in the right-hand side of ( |7.27| l. 
Indeed, we have the following. 

Proposition 7.6. For every A e]0, 1], G(X, ■) e E 1,2 . Moreover, the following state- 



ments hold: (i) E 



J G 2 (X,s)ds 



< d(l + |.T | 2 ), (ii) E 



/ {D Sl G{\s)) 2 dsids 
jo,tp 



< 



C 2 (l + \x \ s ), where C\ and C 2 depend on T, ||cr'|| ^ and ||&'|U, but not on x . 



Proof of Proposition \ 7. 6\ (i) First we set I± = E 



jG 2 (X,s)ds 

.0 



recalling that 



G(M = ^',AZ +3 : )£(A; a ,t) ) ^ ^ =< ^ ^ > 



By Cauchy-Schwarz inequality, we have 

J x < I E J cr 4 (sA 2 , AZ S + a: )ds 

V Lo 
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£ 4 (X;s,t) 



G 



ds 



den 



(7.28) 
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Since a has linear growth, by Lemma [7.2| and using again Cauchy-Schwarz ineguality, 
there is a constant C 1 such that 



h < d (1 + N 2 ) ( E [G den - 8 ] J E [£ 8 (A; s, t)] 



(7.29) 



Conseguently using Proposition 1 7 . 5 1 and Lemma [7~4] the first item of Proposition |7.6| is 
established. 



(ii) We set J 2 = E / (D S1 G(X, s)) 2 dsxds . 
jo.*] 2 

On one hand, by usual Malliavin differentiation rules, we obtain 

D sl G(X,s) = Xa'(sX 2 ,XZ s + x a )a( Sl X\XZ Sl + x )£(X; s u s) £ ^ *' ^ l [0 ,«](ai) 



(7.30) 



+ a(sX ,XZ s +x ) — <j(sX , XZ S + x )£(X- 1 s,t) — — 5 . 

den Gj 



T den 



The right-hand side being measurable with respect to fl x [0,T] 2 , G(A, ■) will belong to 
I 1,2 if (ii) is established. At this point we need to evaluate D Sl £ (A; s,t). From now on, 
for the sake of simplicity, we will only expose the calculations in the case when b = 0. 
In fact, we have 



D Sl £{X;s,t) = £{X;s,t)D Sl (A / a\£X 2 , XZ e + x )dW e - — / (a') 2 (£X 2 , XZ e + x )d£ 



£(X;s,t) Xa'(s 1 X 2 ,XZ Sl + x )t[ a<t ] (s^ 



+ X 2 a( Sl X 2 ,XZ Sl +x ) / t [sut] (e)a"{eX 2 ,XZ e +x )£(X; Sl ,e)dW e (7.31) 



X 3 a( Sl X 2 ,XZ Sl +x ) / l [si ^)(aV0(^A^ + :r o )£(A; Sl ,£)^ 



On the other hand, we get 

t 



D Sl G dm = 2X J aa\iX 2 ,XZ^+x )t [siA {0cr{sA\XZ Sl +x )£{X- Sl ^)£ 2 {X-tt)d^ 
o 

/, 

+ 2 y a 2 ^X 2 ,XZ i+ x Q )£{X;^t)D Sl £{X;^t)di. 



(7.32) 



Therefore, coming back to ( |7.30[ l and using ( |7.31| l, we obtain that 

h < 4 [Ji + J2 + Js] , 



(7.33) 
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where 



J i = E 



J 2 = E 



t s 

ds J ds± 
11 o 



a'(s\ 2 , XZ S + x )a(s 1 X 2 . XZ Sl + x )£(X; s u s) 



£(X;s,t) 



Gden 



ds 1 dsa 2 (sX 2 , XZ S + x )^-^^- — \ (s 1 )a'{s 1 A 2 , XZ Sl + x ) 



G 



[o,t] 2 



den 



+ a( Sl X 2 ,XZ Sl + x ) J l [siit] (e)a"(eX 2 ,XZ e + x )S(X;s u e) 

s 

(j 2 (sA 2 , XZ S + x a ) 



dW e - a' a" (£X 2 , XZ e + x )d£ 



,h = E 



ds^s- 



[o,t] 2 



^ den 



-£ 2 (X;s,t) (D sl G den y 



with D Sl Gden given in l |7.32[ l. In the sequel we will enumerate constants K x to K 20 - all 
those will not depend on x or t, but eventually on T, a and b. We start estimating J\. 
Since a' is bounded, by Cauchy-Schwarz inequality, we have 



J\<KAWa 



t s 

ds J ds 1 a i (s 1 X 2 , XZ Sl + x ) 
o o 



E 



l 8 

J ds J ds 1 £ i (X;s,t) 



o o 



£ 4 (X;s u s) 

^ den 



Since a has linear growth, Lemma [772] and a further use of Cauchy-Schwarz inequality 
imply that, J\ is bounded by 



K 2 (l + \x \ 2 ) (E [G^J) 5 E 



t s 

ds J d Sl £ le (X;s,t) 
n o 



E 



t s 

ds J d Sl £ 16 (X; Sl ,s) 
I ' o 



Therefore, by Proposition [775] and Lemma [774[ we obtain that 

Ji < K 3 (1 + \x \ 2 ) . 
We go on with the analysis of J 2 . Since a 1 , a" are bounded, we have 



(7.34) 



J 2 < KJE 



dsds^isX 2 , XZ S + x )— -4^— — ( 1 + er 2 (siA 2 , XZ Sl + x )M 2 (s, s x \t) 



[0,t] = 



where M(s,s 1 ;t) = J cr"(lX 2 ,XZe + x )£ (A; s ± , tjdWt, t>sVs u is a martingale having all 

s\/ s\ 

moments because of Lemma |7.4| Since a has linear growth, using Cauchy-Schwarz 
inequality and Lemma [772] we get 



J 2 <A' 5 (l + |x | 4 ) E 



ds£ 16 (X]s,i) 



dsds 1 E[M s (s,s 1 ;t)] 



\o,t] 



Then Burkholder-Davis-Gundy inequality Lemma [774] and Proposition |7. 5| imply 

J2<K 6 (l + \x \ i ). (7.35) 
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Finally we treat J 3 . Applying Cauchy-Schwarz inequality, we have 

<j 4 (s\ 2 ,\Z s +x ) c4 



■h < ME 



dsds x - 



G 



'-£\X;s,t) 



den 



E 



y dsds^Ds^den)" 

[o,t] 2 



Since er has linear growth, again by Cauchy-Schwarz inequality and Lemma [772] we get 

J3<K 7 (l + \x \ 2 ) I E[G de 3 I 2 jy"d s E[£ 16 (A; S ,t)] j IE J d Sl (D Sl G den ] 



by Lemma [T4| and Proposition [T3] it follows 



^3 < i^s(l + ko| 2 ) E 



ds 1 (D Sl G d enT 



(7.36) 



Since a 1 is bounded, l |7.32[ l and Jensen's inequality give 



E 



t 

J d Sl (D Sl G den ) 



< K 9 {A 1 + A 2 ) . 



(7.37) 



where 



A x =E 



A 2 = E 



d Sl / d^ 4 (siA 2 ,AZ Sl + ^(£A 2 ,AZ 5 + :z; )£ 4 (A; Sl ,£)£ 8 (A;£,t) 



.0 si 
t t 



y d Si y ^ 8 (eA 2 ,AZ c+2 ; )£' 4 (A;C,t) (D 3l £(X; £, <)) 4 



o o 



if 10 (l + |x | 8 ) E 



Since a has linear growth, Cauchy-Schwarz inequality and Lemma 7.2 imply that Aj is 
bounded by 

i 

t t 

d Sl J d££ 8 (\; Sl ,0£ 16 (\;Z,t) 

.0 si 

Again, by Cauchy-Schwarz inequality and Lemma [7^4] we obtain 

Ax < K lt (l + \x \ 8 ) . (7.38) 
We proceed estimating A 2 . Using Cauchy-Schwarz inequality A 2 is bounded by 



K 12 E 



y d Sl dta 16 {ZX 2 ,\Z s +x ) E y d Sl d££ 8 (\;^t)(D Sl £(\;U)) 8 
jo,t] 2 J |[o,t] 2 

Since a has linear growth, Cauchy-Schwarz inequality and Lemma [7T2] lead to 



A 2 <K 13 (l + \x \ 8 ) E 



d Sl d££ 16 (\;Z,t) 



[o,t] 2 



E 



d Sl dUD ai £(\;Z,t)y 



[o,t]- 
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Lemma |7.4| implies 



A 2 < K 14 (1 + \x \ H ) E 



J d Sl dH(D Sl £(\;^t)y 
[o,«] 2 



(7.39) 



Since a 1 and a" are bounded, using ( |7.31| ) and Jensen's inequality, it follows that 



E 



J d Sl di{D Sl £{X^,t)Y 
[o,tp 



< K u (i?! + R 2 ) , 



(7.40) 



where 



i?i = E 



R 2 =E 



1 d Sl d££ 16 (A;£,t) 



[o,t] 2 



ds 1 ^ ( r 16 ( Sl A 2 ,AZ Sl +a: )f 18 (A;^i) I M 18 ( Sl) £;t) + / £ 16 (X; s u p)dp 



and M(s u £;t) = J a"(p\ 2 ,\Z p + x )£{\;s 1 ,p)dW p , t>&s u is a square integrable mar- 

ev Sl 



tingale, taking into account Lemma [7^4} Again by Lemma [7~4] Ri is uniformly bounded 
in t and x a . On the other hand, using Cauchy-Schwarz and Jensen's inequalities, R 2 is 
bounded by 



K 16 E 



ds^a 32 ( Sl A 2 , AZ S1 + a; )£ 32 (A; £, t) 



[o,t]= 



E 


J d Sl dC 




[o,t] 2 \ 



Since c has linear growth, again by Cauchy-Schwarz inequality, Lemma fO] and Lemma 
|7.4| we get 



R 2 <K 17 (l + \x \ 16 )-E 



d Sl d^M 32 ( Sl ,^t) 



[o,*] 2 



|16\ 



(7.41) 



Burkholder-Davis-Gundy inequality gives 

R 2 < K la (1 + \x \ 
Coming back to l |7.39| l, using ( |7.41| l and l |7.4Q| l, we obtain 

A 2 < K 19 (1 + \x \ 12 ) ; (7.42) 
thus, replacing fl7.38| l and ( |7.42[ l in ( |7.37| l and coming back to ( |7.36[ l, imply 



J 3 < K 20 (1 + |ar | 



(7.43) 



Consequently, substituting ( |7.34[ l, ( |7.35[ l and l |7.43| l in l |7.33| l, item (ii) of Proposition [7^6] 
is established. □ 



Returning to the proof of Proposition |5.1| and substituting in ( |7.27[ l the right-hand 
side of the first and the second item of Proposition |7.6| the inequality l |5.2[ l is verified. 
Finally this concludes the proof of Proposition [57T] □ 
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